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Equation satisfied by electron-electron mutual Coulomb repulsion energy density 
functional. 
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It is shown that the electron-electron mutual Coulomb repulsion energy density func- 
\^q ■ tional V? e [p] satisfies the equation 



K [A] - vz [A_J 
d ' V W (A(r, - p; - lW) 

where p x N (r) and p 1 N _ 1 (r) are iV-electron and (N — l)-electron densities determined 
from the same adiabatic scaled external potential of the iV-electron system at coupling 
strength 7. 
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I. INTRODUCTION 



Density Functional Theory (DFT) 1 is one of the most important tools for the calculation 
of electronic structure and structural properties of solids and molecules^. In all practical 
applications approximations are used^. It is therefore important to examine the exact formal 
properties of functionals that enter into the formalism as an aid to understanding the under- 
lying formulation and to develop improved approximations. Here the Coulomb interaction 
energy functional 
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is examined. 
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\l/ 7 7 ) is the A-electron wavefunction that yields the density pjj and minimizes 
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T and V ee are the kinetic energy and electron-electron repulsion operators respectively and 
7 scales the electron-electron interaction energy. In this paper it is shown that 



vi [ P l N ] - vz [ P l_ x ] 
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where p x N (r) and p' y N _ 1 (r) are A^-electron and (A^— l)-electron densities determined from the 
same external adiabatic potential at coupling strength 7 as discussed below. As a corollary 
it will be shown that 
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II. PROOF 

In the adiabatic connection approach 4-7 of the constrained minimization formulation of 
density functional theory 1 -! 3 -^ the Hamiltonian H 1 for a system of A" electrons is given by 
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T N + yV£ + vl ext [p 



(5) 



Atomic units, H — e — m—1 are used throughout. T is the kinetic energy operator, 

i=\ 



rjnN 



and jV ee is a scaled electron-electron interaction, 

w=iibwh-v (7) 
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The the external potential 

N 
*JV,ext \PN] = Yl V ext (M J 1\) , (8) 

is constructed to keep the charge density fixed at p N (r) , the ground state charge density of 
the fully interacting system (7 = 1), for all values of the coupling constant 7. The external 
potential has the forrrA^ 

<L(M ; r) = (1 - 7) v ux ([p N }; r) 

+ vl(\p N ]; r) - v2([ PN ]; r)+^ xt (M ; r), (9) 

where vl xt ( [pu] ', r) = f ex t (r) is the external potential at full coupling strength, 7 = 1, 
and fextdjojv] ; r) is non- interacting Kohn-Sham potential. The exchange plus Hartree 
potenti&]r^£v ux ([p N ];r), is independent of 7, while the correlation potential i>2([pat]; r) 
depends in the scaling parameter 7. 

The chemical potential 

Mtv = £^v Wxt [Pn]) - ^_i (v c 7 xt M) (10) 

depends on the asymptotic decay of the charge density^ - — , and hence is independent of the 
coupling constant 7^^. In Eq. ( TTUj) E 1 N _ 1 is the groundstate energy of the (N — l)-electron 
system with the same single-particle external potential v J xt {[pn] ', r ) as the iV-electron sys- 
tem: 
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Note that by construction of v2 xt {[pn] ', r ) , Eq. © Pat = p]v is independent of 7, but the 
groundstate density of the (TV — l)-electron system p y N _ 1 , is a function of 7. 
The correlation energy F 7 [pjv-i] is defined as 9 - 
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where 
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Pn-i 



is the Kohn-Sham (TV — 1) independent particle groundstate wavefunction 



that yields the same density as the interacting (TV — l)-electron system at coupling strength 
7. The correlation part of the kinetic energy is given by 
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The virial theorem can be written as 9 



2T 7 [p] + 7 F ee + J rf 3 rV ([p] ; r') (3p (r') + r'.V'p (r')) = 



or as 
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(14) 



T< [p] + F 7 [p] + J rf 3 rV ([p] ; r') (3p (r') + r'.V'p (r')) = (15) 

where the energy functional F 7 [p] is given byii^ 



F 7 [p] = T 7 [p] + 7 V7e [P] 

= r [p] + 7 (F x [p] + f/[p]) + F 7 [ P ] 



(16) 



T 7 [p] and F 7 [p] are the kinetic and correlation energies for the interacting system at coupling 
strength 7 and density p while E x [p] and U [p] are the exchange and Hartree energies at p. 
Take the functional derivative of Eq. ( Tl5l) with respect to i> 7 (r") = i> 7 xt ([p] ; r") and use 
the chain rule for functional derivatives (here and in the following the assumption is made 
that all the functional derivatives are well defined^): 
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The subscript N in 
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indicates that the functional derivative is taken at constant 
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particle number. A useful expression is the Berkowitz-Parr equation^^ which states that 
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and p is the chemical potential, independent of y^^. In addition 



S 7 (r) = S yy (r) 



where the Fukui function^^ 
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where rp is a constant. 

Multiply (II7D by s f^l p] ( ) and integrate over r". With (JU]) and (jTSJ) the result is 
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Since s -, )s f ,-, is symmetric in r and r', Eq. (1261) . with the relations (1231) and (1251) becomes 
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it follows that if p x (r) = A 3 p (Ar) , the uniformly scaled density, 
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Now consider the Schrodinger equation 
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from which it follows that 
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The Levy constrained minimization approach 3 implies that \I/a ({Ai^}) yields A 3 p(Ar) and 



minimizes ( \l/ 
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where, for an iV-electron system 
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with E^ N ( \ 2 v* [pn] ) and Ej f _ 1 ( X 2 v>- [pw] ) the groundstate energies of the N and (N — 1) 

7 

electron systems with the same external potential X 2 v * \pn\ • From Eq. (l3Tj) it follows that 

p\ x = El (A 2 ^ [ Pn ]) - El_, (a 2 ^ [ Pn \) 
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where the last step follows from Eq. ( TTOj) . The chemical potential p^ is independent of ^ ! ' ' 



Combining ([9]), ( )29l) . ( |32l) and ( |34|) yields (suppressing the subscript N for convenience) 
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= 2p - 2i/" ([p] ; r) - r W ([p] ; r) - ^ ([p] ; r) + 7 ^ ([p] ; r) + 7 jU 7 ([p] ; r) ) . (35) 
Upon substitution of ( 135]) into ([2 
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With the aid of (|32|), (fTBD and (|9|) 
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where 

27 [p] = T> [p] - T° [p] . (38) 

From the definition of E% [p] and T7 [p^~-^, Eqs. ([12]) and ([13j) it follows that 

f^ + 7^?([p]; r W([p]; r ) = o, 

and therefore 

2fj,-2 f d 3 rV ([p] ; r') f (r') - / dVr'.W ([p] ; r') / 7 (r') 
= /^'7^r(r') (39) 



The Fukui function is simply the difference between two densities^: 

f N (O = P 1 N (r) - P l_, (r) . (40) 

Note that by construction only the (N — l)-electron density depends on 7. The external 
potential for the N and (N — 1) systems are the same, and hence the virial theorem^ and 
Eq. flUD imply that, 



dVr'.VV^];!-')/^!-') 
= 2T< [p]J + jVZ [pj,] - TUT ^_J - tKI K_J • (41) 

Combining ( 139]) . (J4"T|) and (132]) leads to the main result of this paper: 

VZ [p]v] - Kl K-i] 

= /^^| ! ^ 1 (p)v(r)-p^ 1 (r)). (42) 

III. COROLLARY: RECURSION RELATIONS 

Recursion relations can be derived at coupling strength interaction 7, when the single 
particle external potential is kept fixed, i.e. for an M-electron system 
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5> 7 xt(K];r,). (43) 
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and pit = p]^ (Ug xt [pjv]) ^ s the M-electron density constructed from a groundstate of the 
M-electron Hamiltonian Hj d = T M + jVj}? + iP M ext \pn\ ■ The potential v\ ext [p/v] for a real 
system is the interaction potential between electrons and nuclei. From Eq. ( l4"2j) . 

Kl K] - vz [pl^] 

= /^^| 1 g 1 (p]v(r)-p],_ 1 (r)) (44) 



and 
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Continuing this pattern leads to 

vz \pI\ - vz [pI-m] 

and since V^ [Pi] = 0, it follows that, 

Kl [Plr] = E / ^ ^ [P TJ W-x (0 " ^-x- £ (r)) • (47) 

L=0 ^ °I J N-L \ r ) 

IV. DISCUSSION AND SUMMARY 

For 7 = 0, 1/° e [p] = E x [p] + [/ [pp^. From Eq. (J42]) then follows the equation 
E x [p l N ] + U [&} - E x [ P %_ ± ] - U [^_J = J d s r m^±pM {p]v (r) _ ^ (r) ) 
or 

-^x [pjv] _ E x [P°N-l] 

= [d 3 r S -%^(ph(r)-p% _ l( r)) 

+i | d 3 rrfV (p]v (r) - p ^ (r)) j-^ (p' N (r') - p ^ (r')) . (48) 

This expression has already been derived by Levy and Gorling^. The relation derived here, 
Eq. fj42l) . is generalization of their work, valid at all coupling strengths. Note that in the 
derivation of Eq. ([3]) use is made of Eq. (TTTT) . ( 130|) and (I3TI) . This implies that derivation 
is correct for v-representable densities^^ and the question whether it is correct for general 
densities remains. 

In summary, a relationship between the Coulomb interaction functionals of a many elec- 
tron system at charge densities that differ by one electron, evaluated at the same external 
potential, was derived. The derivation was done for integer numbers of particles, but can 
be extended to fractional particle numbers. As a corollary, it was shown that the Coulomb 
interaction functional can be expressed as a sum over integrals of functional derivatives of 
the Coulomb interaction functionals and charge densities for all densities that differ form 
the total density by an integer. These relations place stringent constraints on the energy 
functionals that appear in density functional theory and it will be difficult for approximate 
functionals to satisfy these equations. 



REFERENCES 

1 P. Hohenberg and W. Kohn, Phys. Rev. B 136, 864 (1964). 

2 J. P. P. Adrienn Ruzsinszky, Computational and Theoretical Chemistry 963, 2 (2011). 

3 M. Levy, Natl. Acad. Sci. USA 76, 6062 (1979). 

4 J. Harris and R. O. Jones, J. Phys. F 4, 1174 (1974). 

5 D. C. Langreth and J. P. Perdew, Solid State Comm. 17, 1425 (1975). 

6 D. C. Langreth and J. P. Perdew, Phys. Rev. B 15, 2884 (1977). 

7 0. Gunnarson and B. I. Lundqvist, Phys. Rev. B 13, 4274 (1976). 

8 W. Kohn and L. J. Sham, Phys. Rev. A 140, 1133 (1965). 

9 M. Levy and J. P. Perdew, Phys. Rev. A 32, 2010 (1985). 
10 A. G6rling and M. Levy, Phys. Rev. B 47, 13105 (1993). 
n R. G. Parr and W. Yang, Density Functional Theory of Atoms and Molecules (Oxford 

University Press, New York, 1989). 
12 R. M. Dreizler and E. K. U. Gross, Density Functional Theory (Springer- Verlag, Berlin, 

1990). 
13 R. Jones and O. Gunnarsson, Rew. Mod. Phys. 61, 689 (1989). 
14 M. Levy and A. GSrling, Phys. Rev. A 53, 3140 (1996). 
15 M. Levy and A. GSrling, Phys. Rev. B 53, 969 (1996). 

16 J. P. Perdew, R. R. Parr, M. L. Levy, and J. R. Balduz, Phys. Rev. Lett. 49, 1691 (1982). 
17 M. Berkowitz and R. G. Parr, J. Chem. Phys. 88, 2554 (1988). 
18 C. C. P. Fuentalba, E. Chamorro, Int. J. Quantum Chem. 107, 37 (2007). 
19 M. Levy and A. Gorling, Int. J. Quantum Chem. 56, 385 (1995). 



10 



